function X=standardize(X)
[~,K]=size(X);
for k = 1:K
    Xk = X(:,k);
    X(:,k) = (Xk - nanmean(Xk))/nanstd(Xk);
end
return